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1 Introduction 



The construction of radioactive ion beams (RIBs) facilities around the world has generated 
a lot of interest in isospin physics [1-3]. These facilities provide the opportunities to 
study the nuclear reactions involving nuclei with neutron or proton excess. These studies 
are helpful in investigating the structure of rare isotopes and the properties of isospin 
asymmetric nuclear matter. The ultimate goal of isospin physics by heavy-ion collisions of 
neutron-rich radioactive nuclei is to explore the isospin dependence of in-medium nuclear 
effective interactions and the equation of state of asymmetric nuclear matter. 

During the last few decades there have been significant activities in exploring the 
isospin effects in collective flow [4-6] and multifragmentation [7,8]. Various studies have 
been done in recent past to investigate the isospin effects in collective flow and in its 
disappearance (at a particular incident energy called balance energy E^) [4-6,9,10]. 
The isospin effects in collective flow have been explained in literature as the competi- 
tion among various reaction mechanisms, such as nucleon-nucleon collisions, symmetry 
energy, surface properties and Coulomb force. The relative importance among these re- 
action mechanisms is not yet clear [4]. To shed light on the relative importance of above 
mentioned mechanisms, in Ref. [9, 10] Gautam et al. have studied the isospin effects 
in E ba i and its system size dependence throughout the range of colliding geometry. The 
study pointed towards the dominance of Coulomb potential in isospin effects at all the 
colliding geometries. Moreover, the study also pointed that the effect of symmetry energy 
is uniform throughout the mass range and range of colliding geometries. So to look for 
an observable which could show the sensitivity to symmetry energy, Sood has studied 
the N/Z dependence of E^/ for isotopic series of Ca [11]. The study shows that the N/Z 
dependence of E^ is sensitive to symmetry energy and shows insensitivity towards the 
isospin dependence of nucleon-nucleon cross section and henceforth the study revealed 
that the N/Z dependence of symmetry energy can act as a probe of symmetry energy. 
To see the system size effects on N/Z dependence of E M , followed by this, Gautam and 
Sood carried out the above mentioned study throughout the mass range at semicentral 
colliding geometry of h/h max = 0.2-0.4 [12]. Their study showed that the sensitivity of 
N/Z dependence of E^ is more for lighter systems. In the present paper, we plan to 
extend the above study to the whole range of colliding geometries. For the present study 
we use isospin-dependent quantum molecular dynamics (IQMD) model [13, 14]. 
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2 The Model 



The IQMD model is an extension of the QMD model [14], which treats different charge 
states of nucleons, deltas and pions explicitly, as inherited from the Vlasov-Uehling- 
Uhlenbeck (VUU) model [15]. The IQMD model has been used successfully for the analysis 
of a large number of observables from low to relativistic energies. The isospin degree of 
freedom enters into the calculations via symmetry potential, cross sections and Coulomb 
interaction. 

In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

fi(r,P,t) = ^ ex P(-[ f -^W] 2 ^) x e M-[p-ti(t)} 2 -^) (1) 

Nucleons are initialized in a sphere with radius R = 1.12 A 1 / 3 fm, in accordance with the 

liquid-drop model. Each nucleon occupies a volume of h 3 , so that phase space is uniformly 

filled. The initial momenta are randomly chosen between and Fermi momentum (pf)- 

The nucleons of the target and projectile interact by two- and three-body Skyrme forces, 

Yukawa potential and Coulomb interactions. In addition to the use of explicit charge 

states of all baryons and mesons, a symmetry potential between protons and neutrons 

corresponding to the Bethe-Weizsacker mass formula has been included. The hadrons 

propagate using the Hamilton equations of motion: 

dfl d(H) dpi d{H) 
dt dpi ' dt d?i 

with 
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xfj(r,p,t)dfdf dpdp. (3) 
The baryon potential V* J ', in the above relation, reads as 

yv (j? — r) = V s 3 kyrme + Vy ukawa + V c 3 oul + V^ m 
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Here t6 = 4C with C = 32 MeV and Zi and Zj denote the charges of the ith and jth 
baryon, and and T^j are their respective T3 components (i.e. 1/2 for protons and 
— 1/2 for neutrons). The parameters /i and ti,....,t± are adjusted to the real part of the 
nucleonic optical potential. For the density dependence of the nucleon optical potential, 
standard Skyrme-type parametrization is employed. 

3 Results and discussions 

We simulate the reactions of Ca+Ca, Ni+Ni, Zr+Zr, Sn+Sn, and Xe+Xe with N/Z 
varying from 1.0 to 2.0 in small steps of 0.2. In particular we simulate the reactions of 
40 Ca+ 40 Ca, 44 Ca+ 44 Ca, 48 Ca+ 48 Ca, 52 Ca+ 52 Ca, 56 Ca+ 56 Ca, and 60 Ca+ 60 Ca; 56 Ni+ 56 Ni, 
62 Ni+ 62 Ni, 68 Ni+ 68 Ni, 72 Ni+ 72 Ni, and 78 Ni+ 78 Ni; 81 Zr+ 81 Zr, 88 Zr+ 88 Zr, 96 Zr+ 96 Zr, 

104 Zr+ 104 Zrj and H0 Zr+ 110 Zr . 100g n+ 100 Sll! H2 Sn+ 112 Sn> 120 Sn+ 120g n) 129^+129^ &nd 

140 Sn+ 140 Sn; and 110 Xe+ 110 Xe, 120 Xe+ 120 Xe, 129 Xe+ 129 Xe, 140 Xe+ 140 Xe, and 151 Xe+ 151 Xe 
at h/h max = 0.2 - 0.4, 0.4 - 0.6 and 0.6 - 0.8. We also use a soft equation of state along 
with the standard isospin- and energy-dependent cross section reduced by 20%, i.e. a = 
0.8 (Jnn e - The details about the elastic and inelastic cross sections for proton-proton and 
proton-neutron collisions can be found in [13, 16]. The cross sections for neutron-neutron 
collisions are assumed to be equal to the proton-proton cross sections. The reactions are 
followed till the transverse in-plane saturates. In the present study we use the quantity 
"directed transverse momentum (pf r )" which is defined as [17, 18] 

1 A 

(pfH^E^G/WKW, (5) 
1=1 

where y{i) and p x (i) are, respectively, the rapidity and the momentum of the i th particle. 
The rapidity is defined as 

y(,-) = ih |W±fiW , (6 ) 

2 E(i)-pS) 

where E(i) and p* z (i) are, respectively, the energy and longitudinal momentum of the 
i th particle. In this definition, all the rapidity bins are taken into account. It is worth 
mentioning that the has the same value for all fragments types [5,19-21]. Further 
the apparatus corrections and acceptance do not play any role in calculation of the E(, a / 
[19,21,22]. 
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In fig. 1 we display the N/Z dependence of E ba i for b/b max = 0.2-0.4 (top panel), 
0.4-0.6 (middle panel) and 0.6 -0.8 (bottom panel). From figure, we find that at all the 
colliding geometries E ba i follows a linear behaviour with N/Z. The slopes are 33, 25, 21, 
18, and 15 (at b/b max = 0.2-0.4), 58, 41, 27, 19, and 18 (at b/b max = 0.4-0.6) and 187, 115, 
67, 39, and 36 (at b/b max = 0.6-0.8) for the series of Ca, Ni, Zr, Xe and Sn, respectively. 
From figure, we find that 

(i) the N/Z dependence of E bat is steeper for the lighter systems as compared to the heavier 
systems at all the colliding geometries, 

(ii) for a particular isotopic series, the N/Z dependence of Eu is more at peripheral 
colliding geometry. 

(iii) and the change in slope is more for lighter systems as compared to the heavier systems 
when we move from central to peripheral colliding geometries. From figure, we see that 
for Ca series, slope increases by almost 400% when we move from central to peripheral 
collisions, whereas for Xe series increase in slope is almost 150%. 

In Ref. [11] Sood has shown that N/Z dependence of E ba i is sensitive to symmetry 
energy and is insensitive to the isospin dependence of nn cross section. The decrease in 
Ebai with increase in N/Z ratio is due to the enhanced role of repulsive symmetry energy 
for higher N/Z ratios. To check the role of symmetry energy on the N/Z dependence of 
E ba i, in Ref. [12] Gautam and Sood reduces the strength of potential part of symmetry 
energy to zero and calculate the E ba i- So here also to see the effect of symmetry energy, 
we reduce its strength to zero and calculate E ba i for two extreme systems of Ca+Ca and 
Xe+Xe series. The results are displayed in Fig. 2 (open symbols). Top, middle and 
bottom panels represent the results for b/b max = 0.2-0.4, 0.4-0.6 and 0.6-0.8, respectively. 
From figure we see that E ba i increases for both the masses on reducing the strength of 
symmetry potential whereas the slope of N/Z dependence of E bat decreases drastically for 
both the systems. From figure, we also see that percentage change (Am(%)= m ~ m3 ^ mmo// ) 
in slope for Ca (Xe) series is 67 (93) at b/b max = 0.2-0.4 whereas it is 64 (89) at b/b max 
= 0.4-0.6. 

In fig. 3, we display the percentage difference of E ba i (AE W (%)= bal - — — *100) 
between calculations without symmetry energy and with symmetry energy as a function 
of N/Z for at bjb max = 0.2-0.4 (top panel), 0.4-0.6 (middle panel) and 0.6-0.8 (bottom 
panel). From fig. we find that 

(i) the percentage difference increases with N/Z for both the system masses at all the 
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colliding geometries indicating that the role of symmetry energy increases with N/Z ratio 

(ii) and the increase is more sharp for Ca series as compared to Xe series at semicentral 
and semiperipheral collisions whereas at peripheral collisions, the increase is almost the 
same for both the masses. 

(iii) Also, for a particular mass, AEm(%) rises with impact parameter and increase with 
impact parameter is more for Xe series as compared to the Ca series. 

In figs. 4 and fig. 5 we display the system size dependence of E ba i for various N/Z 
ratios varying from 1.0 to 2.0 for b/b max = 0.4 - 0.6 and b/b max = 0.6 - 0.8, respectively. 
The results are displayed by open squares, triangles, circles, diamonds, pentagons, and 
left triangles, for N/Z ratios of 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0, respectively. From fig. 
4, we see that E ba i follows a power law behavior (oc A r ) with system size. The power 
law parameter r is -0.50±0.02, -0.46±0.02, -0.41±0.02, -0.34±0.03, -0.33±0.03, and - 
0.22±0.03, respectively, for N/Z ratios of 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0. We find that the 
power law parameter goes on decreasing as we are moving towards asymmetric nuclear 
matter (higher N/Z). This is due to the fact that for higher N/Z ratios, the effect of 
symmetry energy is more in lighter masses (as in Ref. [12]) and thus decreasing by 
larger magnitude in lighter masses which results in less slope for higher N/Z ratio. 

Fig. 5 displays the system size effect of E ba i at peripheral colliding geometry of b/b max 
= 0.6 - 0.8. The behavior of E ba i with system size is similar, except that now the value of 
power law parameter is increased. The value of r now reads as -1.03±0.01, -0.87±0.03, 
-0.81±0.02, -0.72±0.02, -0.63±0.03, and -0.59±0.03,respectively, for N/Z ratios of 1.0, 
1.2, 1.4, 1.6, 1.8, and 2.0. The increase in the value of slope parameter is due to the fact 
that E ba i for lighter systems (like Ca+Ca) changes drastically with impact parameter but 
the change of E ba[ with impact parameter in heavier masses is less [23] which leads in the 
increase of slope parameter at peripheral colliding geometries. 

In Fig. 6 we show N/Z dependence of r for the impact parameter bins of b/b max = 0.2 
- 0.4 (black hexagons) (r values taken from Ref. [11]), b/b max = 0.4 - 0.6 (green hexagons) 
and b/b max = 0.6 - 0.8 (orange hexagons). We see that r decreases with increase in N/Z 
ratio for all the colliding geometries and follows a linear behavior with N/Z having slopes 
0.13, 0.27, and 0.43 for b/b max = 0.2 - 0.4, 0.4-0.6, and 0.6 - 0.8, respectively. We also find 
that change in slope is more when we move from semi peripheral to peripheral geometry 
as compared to that when we move from semi central to semi peripheral collisions. 
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4 Summary 



We have studied the role of colliding geometry on N/Z dependence of balance energy 
(Ey) for isotopic series throughout the mass range. We found that dependence of on 
N/Z ratio is much stronger for peripheral collisions as compared to the central collisions. 
Our study also pointed out that lighter systems showed more sensitivity to the colliding 
geometry towards N/Z dependence of E M . We have also studied the mass dependence of 
E ba i for the N/Z range from 1.0-2.0 for the whole range of colliding geometry. We found 
that the mass dependence of E bat varies with the N/Z ratio. 

This work is supported by a grant from Centre of Scientific and Industrial Research 
(CSIR), Government of India and Indo- French center vide project no-4101-A, New Delhi, 
India. 
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Figure 1: N/Z dependence of for various systems. Various symbols are explained in 
the text. Lines are linear fit. Top, middle and bottom panels represent the results for 
h/b max = 0.2 - 0.4, 0.4 - 0.6 and 0.6 - 0.8, respectively. 
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Figure 2: N/Z dependence of E ba[ for Ca+Ca and Xe+Xe with E sym on and off. Various 
symbols are explained in text. Top, middle and bottom panels represent the results for 
b/b max = 0.2 - 0.4, 0.4 - 0.6 and 0.6 - 0.8, respectively. 
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Figure 3: N/Z dependence of AE W (%) for Ca and Xe series. Various symbols are ex- 
plained in text. Lines are guide to the eye. H 
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Figure 4: System size dependence of for various N/Z ratios at b/b r 
Lines are of power law nature (oc A r ). 
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Figure 5: Same as fig. 4 but for b/b max = 0.6-0.8. 
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